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Abstract 

In this paper we consider a generalization of the classical fc-center problem with capacities. Our goal 
is to select k centers in a graph, and assign each node to a nearby center, so that we respect the capacity 
constraints on centers. The objective is to minimize the maximum distance a node has to travel to get 
to its assigned center. This problem is NP-haid, even when centers have no capacity restrictions and 
optimal factor 2 approximation algorithms are known. With capacities, when all centers have identical 
capacities, a 6 approximation is known with no better lower bounds than for the infinite capacity version. 

While many generalizations and variations of this problem have been studied extensively, no progress 
was made on the capacitated version for a general capacity function. We develop the first constant factor 
approximation algorithm for this problem. Our algorithm uses an LP rounding approach to solve this 
problem, and works for the case of non-uniform hard capacities, when multiple copies of a node may 
not be chosen and can be extended to the case when there is a hard bound on the number of copies of a 
node that may be selected. In addition we estabhsh a lower bound on the integrality gap of 7(5) for non- 
uniform (uniform) hard capacities. In addition we prove that if there is a (3 — e)-factor approximation 
for this problem then P = NP. 

Finally, for non-uniform soft capacities we present a much simpler 11-approximation algorithm, 
which we find as one more evidence that hard capacities are much harder to deal with. 
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1 Introduction 



The fc-center problem is a classical facility location problem and is defined as follows: given an edge- 
weighted graph G = {V, E) find a subset 5 C y of size at most k such that each vertex in V is "close" 
to some vertex in S. More formally, once we choose S the objective function is max^gy miiitig^ d{u, v), 
where d is the distance function (a metric). The problem is known to be NP-hard [il3il . Approximation algo- 
rithms for the A;-center problem have been well studied and are known to be optimal lfT4l[r7l[T8l[T9.| . In this 
paper we consider the /c-center problem with non-uniform capacities. We have a capacity function L defined 
for each vertex, hence L{u) denotes the capacity of vertex u. The goal is to identify a set S of at most k cen- 
ters, as well as an assignment of vertices to "nearby" centers. No more than L{u) vertices may be assigned to 
a chosen center at vertex u. Under these constraints we wish to minimize the maximum distance between a 
vertex v and its assigned center (t){v). Formally, the cost of a solution S is min5cy,|5|=ifc niax^gy d{v, (f){v)) 
such that \{v \ (j){v) = u}\ < L{u) £ S where <j) iV ^ S. 

For the special case when all the capacities are identical, a 6 approximation was developed by Khuller 
and Sussmann |22| improving the previous bound of 10 by Bar-Ilan, Kortsarz and Peleg |2|. In the special 
case when multiple copies of the same vertex may be chosen, the approximation factor was improved to 
5. No improvements have been obtained on these results in the last 15 years. The assumption that the 
capacities are identical is crucial for both these approaches as it allows one to select centers and then "shift" 
to a neighboring vertex. In addition, one can use arguments such as [^] is a lower bound on the optimal 
solution; with non-uniform capacities we cannot use such a bound. This problem has resisted any progress 
at all, and no constant approximation algorithm was developed for the non-uniform capacity version. 

In this work we present the first constant factor approximations for the fc-center problem with arbitrary 
capacities. Moreover, our algorithm satisfies hard capacity constraints and only one copy of any vertex is 
chosen. When multiple copies of a vertex can be chosen then a constant factor approximation is implied by 
our result for the hard capacity version. For convenience, we discuss the algorithm for the case when at most 
one copy of a vertex may be chosen. Our algorithms use a novel LP rounding method to obtain the result. 
In fact this is the first time that LP techniques have been applied for any variation of the fc-center problem. 

While our constants are large, we do show via integrality gap examples that the problem with non- 
uniform capacities is significantly harder than the basic /c-center problem. In addition we establish that 
if there is a (3 — e) -approximation for the A;-center problem with non-uniform capacity constraints then 
P = NP. Such a result is known for the cost A;-center problem [ 8 ] and from that one can infer the result 
for the unit cost capacitated fc-center problem with non-uniform capacities, but our reduction is a direct 
reduction from Exact Cover by 3-Sets and considerably simpler. We would like to note that for the k- 
supplier problem, which can be seen as A;-center with disjoint sets of clients and potential centers, a simple 
proof of (3 — e) approximation hardness under P ^ NP was obtained by Karloff and can be found in 1.18] . 

In all cases of studying covering problems, the hard capacity restriction makes the problems very chal- 
lenging. For example, for the simple capacitated vertex cover problem with soft capacities, a 2 approxi- 
mation can be obtained by a variety of methods |[T5l IT2il - however imposing a hard capacity restriction 
makes the problem as hard as set cover [9]. In the special case of unweighted graphs it was shown that a 3 
approximation is possible H, which was subsequently improved to 2 ifTTTl . 

1.1 Related Facility Location Work 

The facility location problem is a central problem in operations research and computer science and has 
been a testbed for many new algorithmic ideas resulting a number of different approximation algorithms. 
In this problem, given a metric (via a weighted graph G), a set of nodes called clients, and opening costs 
on some nodes called facilities, the goal is to open a subset of facilities such that the sum of their open- 
ing costs and connection costs of clients to their nearest open facilities is minimized. When the facilities 
have capacities, the problem is called the capacitated facility location problem. The first constant-factor 
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approximation algorithm for the (uncapacitated) version of this problem was given by Shmoys, Tardos, and 
Aardal [30 1 and was based on LP rounding and a filtering technique due to Lin and Vitter (26). A long series 
of improvements culminated in a 1.5 approximation due to Byrka |5|. Up to now, the best known approx- 
imation ratio is 1.488, due to Li [25] who uses a randomized selection in Byrka's algorithm [5]. Guha and 
KhuUer |[T6l showed that this problem is hard to approximate within a factor better than 1.463, assuming 
NP ^ DTIME[rf^^°^^''^''\ 

Capacitated facility location has also received a great deal of attention in recent years. Two main vari- 
ants of the problem are soft-capacitated facility location and hard-capacitated facility location: in the latter 
problem, each facility is either opened at some location or not, whereas in the former, one may specify any 
integer number of facilities to be opened at that location. Soft capacities make the problem easier and by 
modifying approximation algorithms for the uncapacitated problems, we can also handle this case |[30ll20l . 
Korupolu, Plaxton, and Rajaraman [23] gave the first constant-factor approximation algorithm that han- 
dles hard capacities, based on a local search procedure, but their approach works only if all capacities are 
equal. Chudak and Williamson [71 improved this performance guarantee to 5.83 for the same uniform ca- 
pacity case. Pal, Tardos, and Wexler [28| gave the first constant performance guarantee for the case of 
non-uniform hard capacities. This was recently improved by Mahdian and Pal [27] and Zhang, Chen, and 
Ye ||3T|| to yield a 5.83-approximation algorithm. All these approaches are based on local search. The 
only LP-relaxation based approach for this problem is due to Levi, Shmoys and Swamy [24] who gave a 
5-approximation algorithm for the special case in which all facility opening costs aie equal (otherwise the 
LP does not have a constant integrality gap). The above approximation algorithms for hard capacities are 
focused on the uniform demand case or the splittable case in which each unit of demand can be served by 
a different facility. Recently, Bateni and Hajiaghayi [4] considered the unsplittable hard-capacitated facility 
location problem when we allow violating facility capacities by a 1 + e factor (otherwise, it is NP-hard to 
obtain any approximation factor) and obtain an O(logn) approximation algorithm for this problem. 

A problem very close to both facility location and A;-center is the k-median problem in which we want 
to open at most k facilities (like in the A;-center problem) and the goal is to minimize the sum of connection 
costs of clients to their nearest open facilities (like facility location). If facilities have capacities the problem 
is called capacitated k-median. The approaches for uncapacitated facility location often work for /c-median. 
In particular, Charikar, Guha, Tardos, and Shmoys ||6l gave the first constant factor approximation for k- 
median based on LP rounding. The best approximation factor for A;-median is 3 + e, for an arbitrary positive 
constant e, via the local search algorithm of Arya et al. [1]. Unfortunately obtaining a constant factor 
approximation algorithm for capacitated /c-median still remains open despite consistent effort. The methods 
used to solve uncapacitated A;-median or even the local search technique for capacitated facility location all 
seem to suffer from serious drawbacks when trying to apply them for capacitated /c-median. For example 
standard LP relaxation is known to have an unbounded integrality gap [6|. The only previous attempts 
with constant approximation factors for this problem violate the capacities within a constant factor for the 
uniform capacity case [6] and the non-uniform capacity case lITOll or exceed the number k of facilities by a 
constant factor [3 |. 

Removing the metric: We employ the standard "thresholding" method used for bottleneck optimization 
problems. We can assume that we guess the optimal solution, since there are polynomially many distinct 
distances between pairs of nodes. Once we guess the distance correctly, we create an unweighted graph 
consisting of those edges uv such that d{u, v) < OPT. We henceforth assume that we are considering the 
problem for an undirected graph G. 

Capacitated /c-Center Problem 

Input: An undirected graph G = {V, E), a capacity function L : F — )• N and an integer k. 

Output: A set 5 C y of size k, and a function : y — )• 5, such that for each u £ S, \(j)~^{u)\ < L{u). 

Goal: Minimize max^jgy distG(f , (t>{v))- 
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By a c-approximation algorithm we denote a polynomial time algorithm, that for an instance for which 
there exists a solution with objective function equal to 1, returns a solution using distances at most c. Note 
that the distance function dist(ii, v), measures the distance in the unweighted undirected graph. 

In the soft-capacitated version S can be a multiset, that is one can open more than one center at a vertex. 
To avoid confusion we call the standard version of the problem hard-capacitated. 

1.2 Our results 

While LP based algorithms have been widely used for uncapacitated facility location problems as well as 
capacitated versions of facility location with soft capacities, these methods are not of much use for problems 
in dealing with hard capacities due to the fact that they usually have an unbounded integrality gap [6il28l. 

For general undirected graphs this is also the case for the capacitated A;-center problem. Consider the LP 
relaxation for the natural IP, which we denote as LPl. We use ?/„ as an indicator variable for open centers. 
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For the sake of presentation we have introduced variables Xu,v for all n, v, even if the distance between 
n and n in G is greater than one. We will use those variables in our rounding algorithm. Furthermore in 
constraints ([T]l and Q we used equality instead of inequality to make our rounding algorithm and lemma 
formulations simpler. In the soft-capacitated version the < 1 part of constraint Q should be removed. 
Note that we are only interested in feasilibity of LPl, and there is no objective function. 

For an undirected graph G = {V, E) and a positive integer 5, by we denote the graph iy, E'), where 
uv G E' iff distG(n, v) < 6. By an integrality gap of LPl we mean the minimum positive integer 6 such that 
if LPl has a feasible solution, then the graph G^ admits a capacitated A;-center solution. As this is usually 
the case for capacitated problems, by a simple example we prove LPl has unbounded integrality gap for 
general graphs. 

Theorem 1. LPl has unbounded integrality gap, even for uniform capacities. 

Proof. Let G' be a graph that consists of two adjacent vertices a, b together with 4 vertices adjacent to both 
of them (see Fig.[T]l. Set uniform capacity L = 4, A; = 3 and consider the graph G which is a disjoint union of 
two copies of G'. Observe that by setting Ua = Ub = 0.75 in each of the copies as well as Xa,v = x^^v = 0.5, 
for all six vertices v, we obtain a feasible solution to the LP relaxation. No matter what 5 we choose, there 
is no capacitated /c-center with L = 4 in the graph G^ . □ 

However, interestingly, if we assume that the given graph is connected, the situation changes dramati- 
cally. Our main result is, that both for hard and soft capacitated version of the A;-center problem, even for 
non-uniform capacities, LPl has constant integrality gap for connected graphs. Moreover by using novel 
techniques we show a corresponding polynomial time rounding algorithm, which consists of several steps, 
described at high level in the following subsection. The actual algorithm is somewhat complex, although it 
can be implemented quite efficiently. 
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Figure 1: A single connected component of the graph G used in the proof of Theorem [T] 

Theorem 2. There is a polynomial time algorithm, which given an instance of the hard-capacitated k- 
center problem for a connected graph, and a fractional feasible solution for LPl, can round it to an integral 
solution that uses non-zero variables for pairs of nodes with distance at most c. 

Corollary 3. The integrality gap of LPl for connected graphs is bounded by a constant, and there is a 
constant factor approximation algorithm for connected graphs. 

To simpUfy the presentation we do not calculate the exact constant proved in the above corollary, but it 
is in the order of hundreds. As a counterposition, for soft capacities in Section[3] we present a much simpler 
11-approximation algorithm, which we find as one more evidence that hard capacities are much harder to 
deal with. 

Theorem 4. For connected graphs there is a polynomial time rounding algorithm, upper bounding the 
integrality gap of LPl by 11 for soft-capacities. 

By using standard techniques one can restrict the capacitated A;-center problem to connected graphs. 

Theorem 5. If there exists a polynomial time c-approximation algorithm for the (soft) capacitated k-center 
problem in connected graphs, then there exists a polynomial time c-approximation algorithm for general 
graphs. 

Proof. Let us assume that we are given a graph G with connected components Ci, . . . , Gt, a capacity func- 
tion L : F— s-N and an integer k. For each connected component Cj using binary search we find the smallest 
value of ki < k for which our black box algorithm finds a solution (if there is no solution for ki = k then 
we set ki = oo). If Yl\=i ki > k then we answer NO, otherwise we return the union of solutions found by 
the black box algorithm. 

To prove correctness of the above algorithm it is enough to observe that if there exists a solution S to 
the A;-center problem in the graph G then for each connected component C, we have ki <\S r\ V{Gi)\ . □ 

Therefore we prove there is a constant factor approximation algorithm for the hard-capacitated A;-center 
problertj^ Our results easily extend to the case when there is an upper bound U{u) of the number of times 
vertex u may be chosen as a center. Constraint [5] should be modified to be < < U{u) to yield a 
relaxation LP2. We can employ the same rounding procedure as discussed for the hard capacity case with 
U{u) = 1. 

The proof of the following theorem is omitted. 

Theorem 6. There is a polynomial time algorithm, which given an instance of the hard-capacitated k- 
center problem for a connected graph, and a fractional feasible solution for LPl, can round it to an integral 
solution that uses non-zero variables for pairs of nodes with distance at most c. 

'with some care, perhaps some of the constants can be improved, however our focus was to show that a constant approximation 
is obtainable using LP rounding. 
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While our constants are large, we do show in Section [5] via integrality gap examples that the problem 
with non-uniform capacities is significantly harder than the basic /c-center problem. 

Theorem 7. For connected graphs the integrality gap ofLPl is at least 5 for uniform-hard-capacities and 
at least 4 for uniform-soft-capacities. 

Moreover in the non-uniform hard-capacitated case, the integrality gap ofLPl for connected graphs is 
at least 7, even if all the non-zero capacities are equal. 

Despite the fact, that the algorithm of f225 for uniform capacities was obtained more than a decade ago, 
no lower bound for the capacity version (neither soft nor hard), better than the trivial 2 — e, derived from the 
uncapacitated version, is known. We beUeve that the integrality gap examples, presented in this paper, are 
of independent interest since they may help in proving a stronger lower bound for the capacitated fc-center 
problem with uniform capacities. 

To make a step in this direction we investigate lower bounds for the non-uniform case. By a reduction 
from the cost fc-center problem [Si one can show that there is no (3 — e)-approximation for the capacitated 
fc-center problem with non-uniform capacities. By a simple reduction from Exact Cover by 3-Sets, in 
Section |6j we prove the same result under the assumption P ^ NP. 

Finally we give evidence that our LP approach might be the proper tool for solving the capacitated 
A;-center problem. The proof of the following theorem shows that when the KhuUer-Sussmann algorithm 
fails to find a solution then in fact there is no feasible LP solution for that guess of distance. The smallest 
radius guess for which the algorithm succeeds, proves an integrality gap on the LP. Considering the result of 
Theorem |7] as we show in Section |4j it follows that for uniform capacities the gap in the analysis is small, 
since our bounds are tight up to an additive +1 error. 

Theorem 8. For connected graphs the integrality gap ofLPl is at most 6 for uniform-hard-capacities and 
at most 5 for uniform-soft-capacities. 

1.3 Our techniques 

We assume that G is connected and that LPl has a feasible solution for the graph G. We call two functions 
X : V X y— >]R_|_ U {0} and y : F— >]R_|_ U {0} an assignment even if {x, y) is potentially infeasible for LPl. 
In other words initially we have a feasible fractional solution, in the end we will obtain a feasible integral 
solution, although during the execution of our rounding algorithm an assignment {x,y) is not required to 
be feasible. Furthermore without loss of generality we assume that for a vertex v with L{v) = we have 
Vv = 0. 

We need to show that there exists a constant 5 such that if for a connected component LPl has a feasible 
solution, then one can (in polynomial time) find an integral feasible solution for G^ . 

Definition 9 ((5-feasible solution). An assignment is called 6-feasible if it is feasible for the graph G^. 

Note that the only difference between LPl's for the graphs G and G^ is constraint (j6|. 

Definition 10 (radius(^ ,^)). For a (5-feasible solution (x, y) to LPl we define a function radius(2; : 
F— )'{0, . . . , (5} which for a vertex u assigns the greatest integer i such that there exists a vertex v with 
distG(f , u) = i and Xu,v > (if no such i exists then radius(^. (u) = 0). 

We give a brief overview of the rounding algorithm described in subsequent subsection of Section [2] 
Initially we start with a 1-feasible (fractional) solution (x, y) to LPl and our goal is to make it integral. We 
perform several steps where in each step we get more structure on the (^-feasible solution but at the same 
time the value of 5 will increase. 
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In Sections 2.ip.4 in four non-trivial steps we round the y-values of a feasible solution. First, in Sec- 



tion 2. 1 we define a caterpillar structure which is a key structure in the rounding process. We show that in 
polynomial time we can find a 5-feasible solution together with a caterpillar structure [P, P') such that all 
vertices outside of the caterpillar structure have integral y-values. In Section 2.2 we define the y-flow and 
chain shifting operations which allow for transferring y-values between distant vertices using intermediate 
vertices on the caterpillar structure. Unfortunately, because the capacities are non-uniform and hard, to find 
a rounding flow for a caterpillar structure we need more assumptions. To overcome this difficulty in the 
most challenging part of the rounding process, that is in Section 2.3 we define a safe caterpillar structure 
and show how to split a given caterpillar structure into a set of safe caterpillar structures (at the cost of 
increasing radius of the (5-feasible solution). In Section 2.4 we design a rounding procedure for a safe cater- 
pillar structure, obtaining a c-feasible solution with integral y-values, for some constant c. We would like to 
note, that for uniform capacities every caterpillar structure is safe, therefore for non-uniform capacities we 
have to design much more involved tools comparing to the previously known uniform capacities case. 

Finally in Section 2.5 we show, that using standard techniques, when we have integral y-values then 
rounding x-values is simple, obtaining a constant factor approximation algorithm. 



2 LP rounding for hard-capacities 
2.1 Group shifting and caterpillar structure 

In the first phase of our procedure we obtain a path-like structure containing all vertices with non-integral 
y- values. We first define the notion of shifting values between variables of LPl relaxation. 

Definition 11 (siiifting). For an assignment {x, y) for the LP, two distinct vertices a,h and a positive 
real a < mm{ya, 1 — yb) such that L(a) < L{b) by shifting a from o to 6 we consider the following 
operation: 

1. Let e = ^; for each v € V let = exa,v, decrease Xa,v by A^, and increase Xh^v by A^,. 

2. Increase y^by a, and decrease by a. 

Lemma 12. Let (x, y) be a 6-feasible solution to LP. Let {x', y') be a result of shifting a from a to b, for 
some a, a, b such that L(a) < L{b), < a < min(ya, 1 — y;,). Then (x', y') is a [5 + distG(a, b))-feasible 
solution and for each vertex v ^ b we have radms(^x' ,y')iv) < i'a'dius(2. (f ) whereas radius(a./ < 
max(radius(^ J,) (a) + distG(a, b), radius (^x,y)it>)) ■ 

Proof. First we prove that {x' ,y') is a ((5 + distG'(a, 6))-feasible solution. Since the sum of all y-values in 
(x, y) is equal to the sum of all y- values in {x', y'), constraint ([T]) of LPl is satisfied. To prove (H of LPl, 
it is enough to consider the variables x'^j x'^, ^ for each vertex v, that is: 

x'b,v = Xb,v + exa,v < yb + ^Va = yb + a = y[ 

For ([3]) of LPl, we only verify v = a and v = b since for other vertices the sum did not change. 

X] = '^i^b,v + (^Xa,v) = Xb^v) + Xa,v) < L{b)yb + eL{a)ya < L{b)yb + eL{b)ya 

= L{b){yb + eya) = L{b)yi 
^'<^,'" = "^i^^,^ - ^^i.f) = (1 - ^) X] - (1 ~ ^)L{a)ya = L{a)y'a 

v&V v&V v&V 
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For each vertex u the sum X^^eF ^v,u is equal to Ylvev ^v,u hence constraint Q is satisfied. Constraints 
([Sjl, (j7]l may be checked directly, since a < min(ya, 1 — yb). 

Since for each vertex u / 6 a variable Xy^u can only be decreased (when v = a), therefore Tadms(^x',y') (^) < 
radius(a, For v = b the new radius may increase but it will not exceed radius(j. (a) + distG(a, b) 

since if x'j^^ > then either > or Xa^u > 0. Consequently constraint (6jl is satisfied for G^' where 
5' = 6 + distG(a, b) and the lemma follows. □ 

Definition 13 (group sliifting). For a (^-feasible solution (x, y) and a set Vb ^ by a group shifting we 
denote the following operation. Assume Vq = {vi, . . . ,V£}, where L{vi) < L{vi-\-i) for I < i < £. As 
long as there are at least two vertices in Vq with fractional y-values, let a be the smallest, and b the greatest 
integer such that Va, Vh G Vq are vertices with fractional y-values. Shift min{ya, 1 — yb) from a to b. 

Lemma 14. Let {x, y) be a 5-feasible solution, Vq be a subset ofV and d = max^^bgyy distG(a, b). After 
group shifting on Vq we obtain a [5 + d)-feasible solution (x' , y'), where there is at most one vertex in Vq 
with fractional y-value and moreover for v £ V \ Vowe have radius(j./ < radius (^x,y){v). 

To make a graph Hamiltonian we use the following lemma known from 1960 ||2T1|291 . 

Lemma 15. For any undirected connected graph G there always exists a Hamiltonian path in G^ and one 
can find it in polynomial time. 

We define a caterpillar structure which is one of the key ingredients of our rounding process. Intuitively 
we want to define an auxiliary path-like tree, where adjacent vertices are close in the original graph G, 
vertices with fractional y- values are leaves of the tree, and all non-leaf vertices have y- values equal to 1. 

Definition 16 (caterpillar structure). By a (^-caterpillar structure for an assignment (x, y) we denote a 
sequence of distinct vertices P = {vi,. . . ,Vp) together with a sequence P' = (vq, . . . , Wp+i) where: 

1. for each i = 1, . . . ,p we have y^, = 1, 

2. for each z = 1, . . . ,p — 1 we have distcCvi, Wj+i) < 

3. for each i = . . . ,p + 1 either v[ = nil orv[£V\ {vj : j = 1, . . . , p}, 

4. for each z = 1, . . . ,p if t;^ 7^ nil then L{vi) > L{v[), < y^/ < 1, distG(wi, v'^) < 6, 

5. if Vq / nil then distG(^^0' ^1) ^ ^ < Vv'^ < !> 

6. if v'p_^_l / nil then distG(t'p+i, fp) < 5,0 < y^>^^^ < 1, 

7. for each < i < j < p + 1 if f ■ / nil and Vj / nil then f ■ / Vj, 

We sometimes omit 6 and simply write "caterpillar structure" when the value of 6 is irrelevant. 

Lemma 17. For a given feasible LP solution (x, y) we can find a 5-feasible solution (x', y') together with 
a 21-caterpillar structure (P, P') such that each vertex v £ V \ (y{P) U V{P')) has an integral y-value in 
{x\ y'), and the first and last element of the sequence P' equals nil. 

Proof. Consider the following algorithm for constructing sets S, S' and a function $ : V^S' . The set S 
will be an inclusionwise maximal independent set in G^ and moreover we ensure that L{<^{v)) > L{v), for 
any v ^V. 
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Figure 2: Example of a (5-caterpillar structure {{vi,V2,vs,V4), {v'q, v[,ml, v'^,v'^, nil)). Vertices connected 
by edges are within distance 6 in the graph G. Note that the sum of y- values over all vertices is integral. 



1. SetFo := l^and^ := S' := 0. 

2. As long as Vo 7^ let w be a highest capacity vertex in Vq- 

• Let f{v) be a highest capacity vertex in Nc[v] (potentially f{v) ^ Vq). 

• Add f{v) to S' and for each u G Ng[Ng[v]] n Vq set $(n) = f{v). 

• Add u to 5 and set Vq := Vq \ Ng[Ng[v]]. 

Observe that each time we remove from the set Vq all vertices that are within distance two from v, hence 
the set S is an inclusion maximal independent set in G^. For this reason vertices in the set S have disjoint 
neighborhoods and moreover by constraints Q and ([2]) of the LPl we infer that for each v ^ V we have: 



1 



(8) 



u&N[v] ueN[v] 

We perform shifting operations to make sure all vertices in the set S' have y- value equal to one. Consider a 
vertex v £ S and the corresponding vertex f{v) chosen by the algorithm. As long as < 1 take any u G 
N[v], u 7^ f{v) such that t/u > and shift mm{yu, 1 — yf{v)) from u to f{v). Note that L{u) < L{f{v)) by 
the definition of f{v) and for this reason shifting is possible. By Lemma [T2| after all the shifting operations 
we have a 3-feasible solution {x,y), since before a shift from u to f{v) we have radius(^ y)(n) < 1, 
radius(a;,j/)(/('f^)) ^ 3 and distG'(?x, f{v)) < 2. Moreover by Inequality ^ we infer, that all the vertices 
in the set S' have y-value equal to one, since otherwise a shifting operation from some u G N[v] to f{v) 
would be possible. 

Observe that by the maximality of the independent set S in G^ the graph G^ [S] is connected, otherwise 
we could add a vertex to S still obtaining an independent set in G^. Moreover for any two adjacent vertices 
u, V G S" in G^[S'], the vertices f{u), f{v) are adjacent in G^[S"]. By the connectivity of G^[S'], the graph 
G''[S"] is also connected. By Lemma 15 we can in polynomial time order the vertices of S' to obtain a 



Hamiltonian path P in G^^ [S']. 

Currently for each vertex v from the set V \ S' we have radius(a; (?;) < 1. For each ?; G S" we use 
group shifting on the set ^^^{f{v)) \ S' . Since 

max distcfa, b) < max distcfa, v) + distcfi', b) < A, 

a,be*-i(/(''))\S" a,6e<I>-i(/(f))\5" 



by Lemma 14 we obtain a 5-feasible solution {x, y) such that all vertices in the set S' have y-value equal 
to one and moreover for each f{v) G S' the set <I>~^(/(f )) \ S' contains at most one vertex with fractional 
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y- value. Let us assume that the already constructed path P is of the form P = {vi, . . . , Vp). We construct a 
sequence P' = (nil, v[, . . . , v'p, nil) where as f • we take the only vertex from ^^^{vi)\S' that has fractional 
y- value, or we set := nilif $~^(tij)\S" has no vertices with fractional y-value. Note that since the way we 
select vertices to the sets S, S' is capacity driven (recall as v we select the highest capacity vertex in Vo and 
as f{v) we select a highest capacity vertex in N[v]), for each vertex u € <I>^^(fj) we have L{u) < L{vi). 
In this way we have constructed a 5-feasible solution (x, y) together with a desired 21 -caterpillar structure 

(P,P')- □ 



As the reader might notice in the above proof we always construct a caterpillar structure with v'^ 







v: 



p+i 



nil. The reason why the definition of a caterpillar structure allows for v'q and v'^^i have non-nil 



values is that in Section 2.3 we will split a caterpillar structure into two smaller pieces and in order to have 



those pieces satisfy Definition 
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we need v'q and v'^j^-^. 



2.2 y-flow and chain shifting 

In the previous section we defined a group shifting operation. Unfortunately we can only perform such an 
operation if vertices are close. In this section we define notions of y-flow and chain shifting which allow us 
to transfer y-value between distant vertices. We will use those tools in Sections [23] and IZ4| 



Definition 18 (y-fiow). For a given assignment (x, y) let S* C F and T (IV he two disjoint sets and let 3" 
be a set containing sequences of the form (a, . . . , vt) representing paths, where a is a positive real, each 

E y is a vertex (for i = 1, . . . ,t), vi & S, vt & T, L{vi) < L{vt) and for i = 2, . . . , t — 1 we have 
Vi ^ S U T, y„. = 1, L{vi) > L{vi). We call {a,vi, . . . , vt) a path transferring a from vi to vt through 
U2, ■ • • , vt-i- We denote f2, . • • , vt-i as internal vertices of the path (a, ui, . . . , uj). 

The set 3" is a y-flow from 5" to T iff: 

• for each v (i S the sum of values transferred from u in 3" is at most y^, 

• for each f G T the sum of values transferred to u in 3" is at most 1 — y^, 

• for each v \ {SVJT) the sum of values transferred through f in 3" is at most 1. 

For a given y-flow 3" from S" to T we define Gj = {V, A) as an auxiliary directed graph with the 
same vertex set as G, where an arc {u,v) belongs to A iff there is a path in 3" containing u and v as 
consecutive vertices in exactly this order. We call the y-flow "J acyclic iff the directed flow graph is 
acyclic. Furthermore we define a function fy : A— )>(0, 1], which for an arc {u,v) assigns the sum of a 
values in all the paths in 3" that contain u and v as consecutive vertices. Moreover by fly : A— )>]R+ we 
denote a function, which for an arc (n, v) assigns the sum of terms L{s)a over all paths from 3" that start 
with a and s G 5 and contain u, v as consecutive elements. Intuitively by /j((n, v)) we denote the fractional 
number of centers that are transferred from u to v, whereas by fl-j{{u, v)) we denote the fractional number 
of vertices (clients) that were previously covered by u and will be covered by v after the shifting operation 
(see Fig.|3]l. 

Now we show that if we are given an acyclic y-flow 3~ then we can transfer y-values using a chain 
shifting method without increasing the radius of vertices by too much. Formal definitions and lemmas 
follow. 

Definition 19 (cliain sliifting). Let 3" be an acyclic y-flow from 5 to T and let (x, y) be a 5-feasible solution. 
Let = {v., A) be the auxihary acyclic flow graph. 
By chain shifting we denote the following operation: 

• For each u,v (iV , set A-u t, = 0. 
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L{si) = 5 




/^((si,a)) = 0.2 



/^((6,ti)) = 0.6 
flA{b,h)) = 1.2. 



L(ii) =4 



fhiia.b)) = 2.6 



/^((^t2)) = 0.2 
flA(.b,t2))^0.4 



L{a) = 12 



/^((S2,a)) = 0.8 
/'^((S2,a))-1.6 



= 6^ 



/^((6,t3)) = 0.2^ 
flA{bM))^l 



\L{t2) = 20 



L{h) - 



Figure 3: The graph Ggr for an acyclic y-flow J = {(0.2, si, a, 6, ta), (0.6, S2, a, 6, ti), (0.2, S2, «, ^, ^2)} 
from S = {si,S2} to T = {hM^h}, where = 0.4, y^^ = = = I, = 0, = 0.8, = 0.1. 
Note that even though each path in 3" has starting point capacity not greater than its ending point capacity 
the vertex G T is reachable from si G 5 in Gg^ despite the fact that L{si) > L{ti). 



• For each arc (u, a) G A in reverse topological ordering of G3-: 

- For each w G V, let A = Xu,vfh{u, a)/{L{u)yu), set Aa,i, = Aa,^ + A and A^,^ = l^u,v - A. 

• For each u,v , set Xu,v = Xu,v + ^u,v- 

• For each s G 5 decrease by E(s,«)eA ^))- 

• For each t G T increase yt by E(„,t)eA hii'^^t))- 

For a directed graph G = {V, A), for a vertex v, we denote A^*"(u) = {n : (u, G A] and A^°"*(t;) = 
{u : (t;,^) G 

Lemma 20. Let (x' , y') be the result of the chain shifting operation on a 6-feasible solution (x, y) according 
to an acyclic y-flow 3" from S to T. If d is the greatest distance in G between two adjacent vertices in 
Gj, then (x',y') is a (5 + d)-feasible solution, and for each vertex v of indegree zero in Gy, we have 
radius(a,' ,^/)(z;) < radius(^ y)(?j), whereas for other vertices v, we have 

vaA\ns{^x' ,y'){v) < max(radius(^^j^)(v), max (radius(a;^j^)(a) + distG(a, w))) . 

Furthermore for each v \ {S[JT) its y-value is the same in (x, y) and {x' , y'). 
Proof. First observe that for any arc (u, v) in Gg- we have the following inequality 

fly{u,v)/L{u)<h{u,v), (9) 

since on each path in 3" the starting vertex has a smallest capacity. Moreover by the definition of a y-flow, 
through each vertex u G 1^ at most units are transferred, hence by Inequality Q we have: 

fHu,a)/L{u)< Yl Mu,a)<yuand (10) 
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We want to show that (x', y') is a (5 + d) -feasible solution. Observe that by the definition of chain 

shifting T^vev v'v = Y^vav V^' ^^'^^ ^^'"tex n e F we have T^veV ^v,u = 'EveV Therefore 

constraints ([T]) and Q are satisfied. Now we check constraint (|7]). Observe that for any u,v & V we have: 



Xu,v + ^u,v > Xu,v - ^ Xu,vfhiu, a)/{L{u)yu) 



(by (fTOb) > - Xu^vVu/Vu = . 



Next we check constraint (j2]l. Consider any u,v ^ V. Observe that if u G V \ {S U T) then either 
x'^ „ = Xu,v ^ Uu = y'u (which happens if u does not belong to any path m3^ or y'^ = 1 > ^ (when u is 
an internal vertex in a path from 3^. On the other hand for s G 5", v G F we have Xs^v < Us and hence by 
Inequality (|9]): 



Similarly for t G T, v G F we have 

aeNg^^it) aeNg^^it) a(iN}^^{t) 

As for constraint Q for a vertex n G y we have: 

E, _ ^-^ >. - X - flyja, u)xa^v s--^ \ - a)xtt,^ 

veV veV veV a£N}V^^(u) ^ ' v£VaeNg^*{u) ^ 

E/, _ v;^ fh{u,a) sr^ sr^ fly{a,u)xa^^ 
"^"'"^ ^ L{u)yu ' ^ ^ L{a)ya 

vev aeAfg"«(«) i-eVaeAfg" («) ^ 

(WQofLP, and C3„ E fif^)- E E^^ 

= L{u)yu- ^ //j(n,a)+ ^ /^3-(a, ^i)/(-i^(a)ya) ^ a;a,i, 
(by(l3)ofLPl)<L(7x)y„- ^ //y(n,a)+ ^ //j(a,n). 

Since for n G 1/ \ (S" U T) by the definition of y-flow EaeAf°"*fM) = EaeAf'" f^l /^3"(0' ^) 

infer T,vev ^u,v < = L{u)y'^. For s G 5 we have N}V-^{s) = and fh{s,a) = /j(s,a)L(s) 

hence J2^^y x',^^ < L{s){ys - EaeAfg";(s) Ms, a)) = L{s)y',. Similarly for t G T we have N^^{s) = 
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and flj{a,t) > fj{a,t)L{t) hence Y.v(iV^'s,v < Ht){ys + Y^aeN^j {t)MO'^t)) = L{t)y[. Therefore 
constraint Q is satisfied. 

Constraint Q may be checked directly by the definition of y-fiow. To check the radius of each vertex in 
(x', y') (that is to verify {5 + d) -feasibility of constraint ([6jl) we observe that if > then either j;^ „ > 

or Xa,v > for some a G N^q {v). 

□ 



2.3 Separable caterpillar structure 



If we knew that in the caterpillar structure (P, P') produced by Lemma 17 the capacity of each vertex in P 
is not smaller than the capacity of each vertex in P' then we could skip this section. Unfortunately some 
vertices of V{P) may have smaller capacity than some vertices of V{P') and for this reason we define the 
notion of dangerous, safe and separable caterpillar structures. 

Definition 21 (safe, dangerous). For a caterpillar structure y = (P = {vi, . . . , Vp), (wq, . . . , Wp+i)), by 

r(J') C V{P) we denote the set containing all vertices Vi, such that there exist {) < < i < ii < p + 1, 
such that v[^^ / nil, L{v[^) > L{vi) and / nil, L{v[^) > L{vi). 

A caterpillar structure T is safe if r(y) = and dangerous otherwise. 

Definition 22 (separable). Let {x,y) be a (5-feasible solution and let J" = (P = {vi,...,Vp), 
P' = {v'q, . . . , fp_|_i)) be a dangerous caterpillar structure. We call T separable iff there exists I < i < p 
such that Vi G r(J'), L{vi) = min^gp(y) L{v) and either: 



• Si > \S2] — 5*2, where S2 = Yl Uv'. 'Si is the sum of values (l — y^) where v ^V,v 
Vj, L{v) > L{vi) for some i < j < p + \, or, 

• Si> [52] — ^2, where ^2 = X] j=o,...,i-i y^,/ and Si is the sum of values (1 — y^) where v ^V,v 

J 

Vj,L{v) > L{vi) for some < j < i. 



We call such i as above a witness of separability of J". A caterpillar structure that is not separable is called 
non-separable. 

The intuition behind the sums Si, S2 is as follows. The sum 5*2 contains all the y- values of vertices of 
P' to the right (or left) of i. Since we want to round all the y-values of vertices of P', if we want to split the 
caterpillar structure (P, P') by removing the edge VjVj+i (or Vi^iVi), we need to send [^2] — 5*2 units of 
flow to the part that does not contain vi, in order to make the sum of y- values over all the leaves in both new 
caterpillar structures integral. That is to satisfy ([8]) of Definition [16] In Si we sum over all vertices, that can 
potentially receive flow if we start a path at Vi, and the value (1 — y^) is the y- value a vertex v may receive. 

An example of a separable caterpillar structure is depicted in Fig. [4] Observe that a non-separable path 
structure may be dangerous as in Fig. [5] 

Lemma 23. Let J" = ((f 1, . . . , Vp), {v'q, . . . , fp+i)) be a dangerous caterpillar structure and let i be an 
index such that Vi G r(T) and L{vi) = min^gp^y-j L[v). Moreover let j be an index such that Vj 7^ 
ml,L{v'j) > L{vi). Then for any a G [min(i,j),max(z, j)] we have L{va) > L{vi). 

Proof. Assume the contrary, that is assume that such a exists. Clearly a ^ i (since then L{va) = L{vi)) and 
also a ^ j because L{vj) > L{v'j) > L{vi) (where the first inequality follows by (|4|) of Def. 16 1. W.l.o.g. 



assume that i < j (the other case is symmetric). By the definition of F (Def. [2T] ) we infer thatthere exists 
< io < i such that v'-^ ^ nil and P(w^y) > P{vi) > L{va). Consequently Va G r(J') which contradicts 
the fact that Vi is a smallest capacity vertex in r(3'). □ 
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Lemma 24. Let 7 = ((vi, . . . , Wp), (wq, • • • , Vp+i)) be a dangerous non-separable caterpillar structure and 
let i = min^gp(y) L{v). For I = {i : < z < p + 1 A t; • / nil A L(f ■) > £} we have Y^ieii^ ~ y-"0 < ^• 

Proof. Consider any Vi G r(y) such that L{vi) = £. Let /i = / n [0, -i - 1] and /2 = / n [i + l,p + 1] 
(note that / = /i U /2). We know that Vi is not a witness of separability hence each of the two sums 5i in 
Definition [22|is strictly smaller than 1, since otherwise we would have S*! > 1 > [^2] — 82- Consequently 

Eje/i(l - y^^) < 1 and similarly T.iehi'^ - Vv^) < □ 



L{vi) = 5 
• 



L(V2) = 3 

= 1 



Liv-s) = J 

— 



L(v4) = 9 

yi-4 = 1 



L(«5) - 3 

■-->• 



L(«6) = 5 



L{vr) - 2 

= 1 



0.2 

= 10 



0.1 



0.3 



0.3 



0.2 



0.2 



0.1 



0.2 ' 



0.8 



yv, = 



0.9 



0.7 



y„;i = 0.3 
L{v'^) = 2 



Figure 4: A separable caterpillar structure (P = (wi, . . . ,V7),P' = {v'q, w'l, nil, ^2,^3, nil, Ug, Ug, nil, nil)), 
where r((P, P')) = {vi,V2,V5} (note that r((P, P'))> since = nil). By dashed edges an acyclic 
flow 3" = {(0.1, V2, vs, V4, ^4), (0.2, V2, W3, V4, ^5, ^6, v'g)} from {V2} to {w^, v'^} is marked with values 
printed in the middle of each arc. 



In the following lemma we use the procedure SEPARATE(y = (P, P')) which given a 5-caterpillar 
structure produces a set of non-separable (5-caterpillar structures. The pseudocode of SEPARATE is given as 
Algorithm [T] and at high level it performs the following steps: 

1. (Lines 1-2) If (P, P') is non-separable then return T. 



2. (Line 3) Otherwise let i be a witness from Definition 22 with the smallest value of L{vi) and assume 
that i is a witness due to the first conditition in the definition (the other case is symmetric). 

3. (Lines 4-5) If the sum S2 is integral, then run the procedure recursively on two caterpillar structures 
{{vi, ...,Vi), {v'q, . . . , nil)), ((ui+i, . . . , Vp), (nil, . . . , v^+i)) and return the set of obtained 
non-separable (5-caterpillar structures. 

4. (Lines 6-9) Let I = {j : i < j < p + 1, Vj nil, L{vj) > L{vi), that is the set of indices used in the 
definition of the sum 5i and denote / = {ii, . . . , where ii < . . . < ij. and let be the smallest 
integer such that X^jg/ j<^Q (1 — Uv'. ) ^ ['521 — S2 (such Iq exists since i is a witness of separability). 

5. (Lines 10-11) Construct an acyclic y-flow 3" from {vi} to {v'^ '■ j ^ I,j ^ ^o}- That is for each 
j £ I, j < £0 add to 3" a sequence (1 — y^/ , Vi, fj+i, . . . ,Vj, v'-) and for j = Iq add to 3~ a sequence 
(a, Vi^Vi+i, veo,v'^^^), where a = ([^2] - ^2) - Eie/,j«?o(l ~ Vv'^) ^^^^ Fig.fflfor illustration); 
perform chain shifting according to 3~. 

6. (Lines 12-18) If i < p then run the procedure recursively on the caterpillar structure ((I'i+i, . . . , Vp), 
(nil, . . . , ^'p+i))' where Vj = Vj if y„/ < 1 and v" = nil otherwise for i + 1 < j < p + 1 (if 
i = pwe already know that 2/^'^^ = 1)- 
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7. (Lines 19-21) If i = 1 then perform group shifting on the set {v'q, vi,v[} \ {nil} and return. 

8. (Lines 22-23) If v'^ / nil then shift min(y^/ , 1 — y^- ) from y^t to y^^ (by the properties of a caterpillar 
structure we know that L{vi) > L{v[)). 

9. (Lines 24-29) If = 1 then run the procedure recursively on the caterpillar structure {{vi, . . . , Vi), 
{v'q, . . . , v'^_i,u, nil)), where u = v'^if v'^ nil, y^i > 0, or n = nil otherwise. 

10. (Lines 30-31) Otherwise (if y^. < l)runtheprocedurerecursively on the caterpillar structure {{vi, . . . ,Vi-i), 

{v'Q,...,v'i_^,Vi)). 



Lemma 25. For a given feasible LP solution (x, y) we can find a 68-feasible solution (x' , y') together with 
a set of vertex disjoint non-separable 21-caterpillar structures S such that each vertex v outside of the set 
has an integral y-value in {x' , y'). Furthermore for each vertex v that belongs to some caterpillar structure 
from S we have radius(^/y)(f ) < 47. 

Proof. First we use Lemma[T7]to obtain a 5-feasible solution (x', y') together with a 21-caterpillar structure 
(i-Q) ^o) such that each vertex outside {Pq, Pq) has an integral y-value in (x', y'). Next we run the procedure 
Separate on the caterpillar structure (Pqi ^o) using (x', y') as the assignment. 

Let us prove, that given a (5-caterpillar structure the procedure Separate indeed returns a set of non- 



separable (5-caterpillar structures. By Lemma 23 the set 3" created in Line 10 of the procedure Separate 



is indeed an acyclic y-flow. Furthermore caterpillar structures created in Lines 5, 18, 29, 31. are indeed 6- 



caterpillar structures since they satisfy all the conditions of Definition 16 Moreover observe that if a vertex 
V belongs to r(J''), where CP' is a caterpillar structure created in one of Lines 5, 18, 29, 31, then v G r(J'). 

Since some vertices of the caterpillar structure created in Lines 18, 29 and 31 may have increased radius 
we need to argue why the radius does not grow too much in the subsequent recursive calls. 

Let S be the set of non-separable 21-caterpillar structures, which is a result of Separate(Po) ^o)- 
Observe that the caterpillar structures are vertex disjoint and moreover each vertex v which does not belong 
to any caterpillar structure from S has an integral y-value in (x', y'). 



Due to Lemmas 12|14|20 when we modify the assignment (x', y') by shifting, group shifting or chain 



shifting we satisfy all constraints of the LP possibly except To show 68-feasibility of (x', y') we prove 
that for each caterpillar structure ((wi, . . . , Vp), {v'q, . . . , fp_(_i)) used as an argument of the procedure Sep- 
arate there exist two indices < < Jo < P + 1 such that: 

• for each i G [l,p] we have radius(^/y)(fi) < 5 + 21 ■ [i < iq] + 21 ■ [i > jo], 

• for each i e [0,p+ 1] if f • / nil then radius(a./^y)(u^) < 5 + 21 ■ [i < io] + 21 ■ [i > jo], 

• for each i G [0, io) U {jo,p + 1] such that v'- ^ nil for any v G r(CP) we have L{y'^ < L{v). 

In the formulas above we use the Iverson's bracket notation, that is [0] equals 1 when (j) is true, and 
otherwise. We call a caterpillar structure satisfying the above properties a good caterpillar structure. The 
initial caterpillar structure {Pq, Pq) is good since we can set zq = and jo = P + 1- Hence we assume 
that a caterpillar structure {P, P') given as an argument to the SEPARATE procedure is good and we want to 
show that all recursive calls are given good caterpillar structures. If the procedure exits in Line 2 then it is 
non-separable, thus we assume that {P, P') is separable and due to the symmetry w.l.o.g. we may assume 



that i is a witness satisfying the first condition of Definition 22 Clearly caterpillar structures constructed in 
Line 5 are good caterpillar structures. Now we investigate caterpillar structures constructed in Lines 18, 29 
and 31. Observe that due to the definition of a good caterpillar structure we can prove that io < i < jo- 
Indeed if i G [1, io] then there is no < i' < i with v'^ / nil and L{v'j) > L{vi) and similarly for i G [jo,p] 
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Algorithm 1 Separate 



Input: A 5-caterpillar structure (J* = (P, P')). 

Output: A set of non-separable 5-caterpillar structures S, such that each vertex v which belongs to T and 
does not belong to any caterpillar structure of S, has integral y- value. 

1: if (P, P') is non-separable then 

2: return {CP} 

3: i ^ a. witness from Definition 22 with the smallest value of L{vi) 

/*assume that i is a witness due to the first condition*/ 
/*in the definition, since the other case is symmetric*/ 

4: if 5*2 is integral then 

5: return Separate((?;i, . . .,Vi), (w^, . . . ,t;^,nil)) U SEPARATE((z;j+i, . . .,Vp), {ml,Vi^^, . . ■,Vp^^)) 
6: S ^ 

7: I ^{j :i<j <p+l,v'j^ nil, L{v'j) > L{vi)} 

/*I is the set of indices used in the definition of the sum Si */ 
8: Denote I = {ii, . . . , ir}, where ii < . . . < ir 

9: io ^ the smallest integer such that J2jei j<eo^^ ~ ^ — 1^*^21 — •5*2 

/*such ^0 exists since i is a witness of separability*/ 
10: Construct an acyclic y-flow 3" from {vi} to {vj j ^ I,j < ^o}- That is for each j ^ I, j < £q add to 3" 
a sequence (1 — y^'_,Vi, I'j+i, . . . ,Vj,Vj) and for j = Iq add to 3~ a sequence (a, Vi, ViJ^i, . . . , u^p, u^^), 
where a = {\S2] - S2) - Y.jei,j<eo^^ ~ Vv'^) (see Fig. |4]for illustration). 
11: Perform chain shifting according to 3~. 



12: if i < p then 



13: 
14: 
15: 
16: 
17: 
18: 
19: 
20: 
21: 
22: 
23: 



for j=i-i- 1 to p-i- 1 do 

if / nil and y^i, < 1 then 



/*ifi = pwe already know that Vv' = 1*/ 



else 



Vj ^ nil 



S = S U SEPARATE((wi+i, . . . , Vp), (nil, . . . , Vp+l))' 

if i = 1 then 

Perform group shifting on {7;q, vi,v[} \ {nil} 
return S 
if 7^ nil then 

Shift min(y„/, I - VvJ from y„/ to y^^ 

/*by the properties of a caterpillar structure we know that L{vi) > L{v'-)*/ 

24: if y^. = 1 then 

if v[ / nil and > then 



25: 
26: 
27: 
28: 
29: 
30: 
31: 
32: 



u = v'i 



else 

u = nil 

S ^ S U Separate((?;i, . . . , Vi), (vq, . . .,v[_]^,u, nil)), 
else 

S ^ S U Separate((ui, . . . , Vi-i), {v'q, v'i_^,Vi)) 
return S 
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there is no i < i' < p + 1 with 7^ nil and L{v'j) > L{vi). Note that for Iq defined in Line 9 we have 
^0 < Jo since due to the last property of good caterpillar structure / C [i + 1, jg]- Furthermore for each 



j G [i + 1,^0) either v'^ = nil or L{v'-) < L{vi) or after chain shifting by Lemma 20 the vertex v'- will 



have integral y-value equal to one and will be not be included in caterpillar structuresmindled recursively. 
Consequently the caterpillar structure constructed in Line 1 8 is a good caterpillar structure. Finally consider 
caterpillar structures created in Lines 29 and 3 1 and observe that a pair of indices io , i satisfies the properties 
of a good caterpillar structure. 

Since all the caterpillar structures from the set S are good caterpillar structures the last part to show 
before we claim that {x',y') is a 68-feasible solution is that vertices from X = {V{Po) U V{Pq)) \ 
(U(pp')eS ^(^) ^(^')) have radius bounded by 68 (by the definition of a good caterpillar structure 
all vertices from U(pp')es ^(^) ^i^') have radius bounded by 47). The only possibilities for a vertex 
V to become a member of X is when its y- value becomes integral, which may happen in Line 17 (when 
y' , = 1), Line 20 (for each non-nil vertex in {v'q, vi, v'l}), Line 23 (when y' , = 0). However in all cases 

because iq < i < jo and by Lemma 20 we infer that radius(3;/ ,^/) (v) < 68. □ 




Figure 5: A dangerous caterpillar structure {P = (f 1, . . . , v^), P' = (?;q, w^, nil, V2,v'^,v'^, v'^, v'q, Vj, nil)), 
where T{{P, P')) = {v-^, v^}. The caterpillar structure is non-separable because both for i = 3 and i = 5 
in Definition [22] the sum Si is at most 0.6, while [^2] — S2 is equal to 0.9. By dashed edges an acyclic 
flow 3~ = {(0.6, vs, V2, vi,Vq), (0.4, V3, V4, v^, vq, vj, Vj)} from {v^^} to {vq,Vy} is marked with values fj- 
printed in the middle of each arc. 



In the following lemma we transform non-separable caterpillar structures into safe caterpillar structures. 

Lemma 26. There exist constants c,6 such that for a given feasible LP solution {x,y) we can find a c- 
feasible solution (x', y') together with a set of vertex disjoint safe 5-caterpillar structures S such that each 
vertex v outside of the set has an integral y-value in (x' , y'). 



Proof. We use Lemma 25 to obtain a set § of vertex disjoint non-separable 21 -caterpillar structures. Our 
goal is to transform each dangerous caterpillar structure in S into a safe caterpillar structure. 

Consider a dangerous non-separable iJo-caterpillar structure J" = ((wi, . . . , Vp), (v'q, . . . , fp+i)) € S and 
let Va be a minimum capacity vertex in r(T). Moreover let / = {i : < i < p + 1 A f • 7^ nil A L{v'j) > 
L{va)}. Construct any acyclic y-flow which sends min(l, J2iei(^ ~ Vv')) from {va\ to {v',- : i ^ 1} (see 



Fig. [5]). Such flow always exists due to Lemma 23 



Let Y = {va, ^a, I'a-i} \ {^^^ perform group shifting on Y (note that a > 1, since Va G r(y)). Re- 
place y in S with the (2(5o)-caterpillar structure {{vi, . . . ,Va-i,Va+i, ■ ■ ■ ,Vp), 
{v'q, . . . , v'^_2, u, v'^_^_i, . . . , Vp_^_i)), where as u we set the only vertex from Y with fractional y-value af- 
ter group shifting or we set u = nil if all vertices in Y have integral y- values. We need to argue, that when 
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It 



u 7^ nil, we have L{u) < L{va~i), in order to satisfy Q of Definition 16 Observe, that if ^ nil, then 
L{v'a) < L{va), and similarly if v'^_i / nil, then L{v'^_^) < L{va-i)- Hence to show L{u) < L{va-i 
is enough to show L{va) < L{va~i), but this follows from Lemma [23| since Va G r(J'). 

Note, that each caterpillar structure will be modified according to the above procedure at most twice, 
since after one iteration the sum ^^^/(l 
we have Y.iaii'^ - Vv^ 

(5-caterpillar structure together with a c-feasible solution. 



2/^' ) either equals zero or decreases by one, and by Lemma 
< 2. Consequently by Lemmas 
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14 



24 



we obtain the desired set of vertex disjomt 

□ 



2.4 Rounding safe caterpillar structures 

In this section we describe how to round the c-feasible solution (x', y') using the set of vertex disjoint safe 



caterpillar structures S from Lemma 26 In order to do that we introduce a notion of rounding flow which is 



a special kind of y-flow defined for a caterpillar structure. 

Definition 27 (rounding flow). For a caterpillar structure (P, P') and an assignment (x, y) we call 9" a 

rounding flow iff 3" is a y-flow from 5 to T where SVJT = V{P'), for each t>- E 5 we have fy{{v[,Vi)) = y^i, 
and for each ■ € T we have /^((f j, u^)) = 1 — y^/. Furthermore each flow path from 3" can not go through 

a vertex from y \ (F(P) U 

In order to obtain a rounding flow for each vertex of V{P') (which by definition have fractional y- 
values), we have to decide whether it will be a source (member of S) or a sink (member of T). After chain 
shifting according to 3~ all sources should have y-value equal to zero whereas all sinks should have y- value 
equal to one and consequently all vertices from the caterpillar structure will have integral y-value. In the 
following lemma we show that for each non-separable caterpillar structure we can always find a rounding 
flow in polynomial time. 

Lemma 28. For any safe 5-caterpillar structure {P, P') and an assignment (x, y) there exists a rounding 
flow 3" such that for any two adjacent vertices in G^r their distance in G is at most 5. Furthermore we can 
flnd such a rounding flow in polynomial time. 

Proof. We present a recursive procedure which constructs a desired rounding flow. Note that some recursive 
calls of the procedure might potentially involve infeasible assignments (x', y'), however we prove that if the 
initial call gives the procedure a safe (5-caterpillar structure, then as a result we obtain a valid rounding flow. 

Let us describe a procedure which is given a caterpillar structure (P, P') together with an assign- 
ment y (the procedure does not need the x part of an assignment). Denote P = {vi, . . . , Vp) and P' = 
{v'q, . . . , v'p^i). If V{P') = then we simply return the empty rounding flow. Otherwise let i be the small- 
est integer such that the sum of y-values of X = {v'q, . . . , v'^} \ nil is at least one (such i always exists 
since the sum of all y-values in V{P') is integral by ([s]) of Def. 16 1. Note that since all vertices in V{P') 



have fractional y-values we have i > 0. Let < io < « be an index such that v'^^^ ^ nil and v'^^^ has the 
biggest capacity in X. Let a = Ylv(^x V-"- If a = 1 then we recursively construct a rounding flow 3" 
from to T for a smaller caterpillar structure ((ui+i, . . . , Vp), (nil, w^^^, . . . , ^^p+i)) and (i) add to S the 
set of vertices X \ {v[^^ (ii) add to T the vertex v[^^ (iii) for each v'^ ^ X \ {v[^^} add to 3" a flow path 
{yy',,Vj,Vj, . . . , Vig,v'^^). In this case we return 3~ as a desired rounding flow for (P, P'). Hence from now 
on we assume a > 1 and q — 1 < y^,/. Consider two cases: io < i and iq = i. 

i 

First let us assume that io < i. We store z := y^' and temporarily set y^i = a — 1. Next recur- 

i i 

sively construct a rounding flow 3" from S C V{P") to T C V{P ) for a smaller caterpillar structure 
{{vi, . . . , Vp), P"), where P" = (nil, v'^, . . . , v'^^i) (note that the sum of y-values in P" is integral). Now 
consider two cases: 
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if f • E 5 then: (i) add to S vertices from X \ {v'^jV^^} (ii) add to T the vertex v^^^ (iii) for each 
VjeX\ {fig, v^} add to 3" a flow path {y^i^.v'^Vj, . . . , t;^^, v^^) (iv) add to 3" a flow path {z — 
y^>,Vi,Vi,..., Vig,Vi^) (v) set y„/ := z (vi) return 3". 

if ■ G r then: (i) add to S vertices from X \ {f^, t'i,,} (ii) add to T the vertex v'-^_^ (iii) out of the flow 
paths in 3~ that end in u leave only that many, that send exactly 1 — z units of flow and reroute the rest 
paths to v'^^ through vertices Vi-i,Vi-2, ■ ■ ■ ,Vig (iv) for each VjeX\ {f , v^} add to 3" a flow path 

^, Jj, . . . ,Vig, Wj^ , 



iVv' ,v'j,Vj,..., Via > ^^io ) '"^^"'■^ ^ 



Now assume that io = i. We create a smaller caterpillar structure {{va,ViJ^i,Vi+2T ■ ■ 
(nil, v'^, v[j^i, . . . , fp+i)), where Va, are two newly created vertices with := a — 1 and L{v'^) := 
L{va) := L{v'-_^), where v'-_^ is the second biggest capacity vertex in the set X. Next run recursively our pro- 
cedure on the newly created caterpillar structure to obtain a rounding flow 3" from S to T. Again, consider 
two cases: 

• ifv'^e S then: (i) set S := {S \ {v'J) U {X \ {v'^}) (ii) set T := T U {v'^} (iii) change in J all the 
paths that start in v'^ to start in X \{vl} (iv) add to 3" paths that start in X and transfer 1 — y^/ units 
of flow from X to v'^ (v) return 3". 

• ifv'^G Tthen: (i)setS' := SU{X\{v'i,v'iJ) (ii)setr := (T\{<})U{u^,t;^J (iii) reroute some of 
the flow paths from 3^ that end in v'^ to that transfer exactly 1 — y^i units of flow to Vii (that is remove 
v'a as the last vertex on those paths and extend the paths by Vi,v'^) (iv) reroute all the remaining flow 
paths in 3" that end in v'^ to v'^^ (that is remove v'^ and extend those paths by Vi, Wj-i , . . . ,Vi-^,v^^) (v) 
for each VjeX\ {v^, } add to 3~ a flow path {y^i ,Vj,Vj, . . . , Vi-^ , ) (v) return 3". 

Finally we prove that if the procedure receives a safe caterpillar structure then it returns a desired round- 
ing flow. The only property of the rounding flow that needs detailed analysis is the assumption that each 
internal vertex of a flow path has capacity not smaller than its the capacity of its starting point. Let us assume 
that there exists a path in 3" that starts in v'^, goes though vi, and ends in v'^, where L{v'^) > L{v'^) > L{vb)- 
This contradicts the assumption that T is safe because G r(J'). □ 



The following theorem summarizes Sections 2.1[ 2.2\ 2.3\ 2.4| 



Theorem 29. For a connected graph G, ifLPl has a feasible solution then we can find a c-feasible solution 
with integral y-values. 



Proof. Using a feasible solution to LPl, by Lemma 26 we obtain a c-feasible solution {x' , y'), together with 



a set of vertex disjoint safe (5-caterpillar structures S, such that vertices that do not belong to any caterpillar 



structure in S have integral y-value in [x' , y'). Next by Lemma 28 for each 5-caterpillar structure (P, P') G S 



we find a rounding flow 3"(p p/). Finally for each J-caterpillar structure (P, P') we perform chain shifting 



with respect to 'Jjppiy and by Lemma 20 we obtain a c'-feasible solution (x", y") to LPl. 



By Lemma 26 vertices outside of Shave integral y-value in (x', y'). Moreover by Definition 
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after 



chain shifting all the vertices in each caterpillar structure of S have integral y-values in (x", y"). □ 
2.5 Rounding x-values 



In this section we show how to extend Theorem 29 to obtain not only integral y-values, but also integral 



x-values. The following lemma is standard (using network flows). 

Lemma 30. Let (x, y) be a 5-feasible solution such that all y-values are integral. There is a polynomial 
time algorithm that creates a 5-feasible solution which has both x- and y-values integral. 



As a consequence of Theorem 29 and the above lemma the proof of Theorem|2]follows. 
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3 Soft capacities 



This section is devoted to the soft capacities variant, where one can open an arbitrary number of centers in 
a node. We present an 11-approximation algorithm for the soft capacitated version of the A:-center problem 
with non-uniform capacities, i.e. we prove Theorem |4] Let us recall that for this problem the LP relaxation 
for the natural IP, which we denote as LPl has the following form. 



E«6vy« = ^; (12) 

Xu,v <yu yu,v £V (13) 

E^ev ^".f ^ L{u)yu Vn G F (14) 

Euev^u,v = i yvev (15) 

Xu,v = Vu, f G V distciujv) > 1 (16) 

Xu,v>0 'iu,veV (17) 

The following lemma is proved by KhuUer and Sussman. 

Lemma 31 ( II22II ). For a connected graph G = {V, E) one can in polynomial time construct an inclusion- 
wise maximal independent set S in G^, such that [S] is connected. 

The proof of the main theorem of this section is inspired by the algorithm of fTI\ for uniform soft 
capacities. 

Theorem 32. There is a polynomial time algorithm, which given an instance of the soft-capacitated k- 
center problem for a connected graph, and a fractional feasible solution for LPl, can round it to an integral 
solution that uses non-zero Xu,v variables for pairs of nodes with distance at most 11. 



Proof. First, we construct the set S using the algorithm from Lemma 31 Observe, that by constraints ( 15 1 
and ( [T3] ) for each v £ V we have 

ueN[v] ueNlv] 



Since S is independent in G^, by constraint ( 12 1 we infer that 15"! < k. 

We prove the theorem in two steps. First, we construct a function (j) : V ^ S, such that for each v £ V 
we have distciv, (p{v)) < 5 and for each s G S" we have |(/)^^(s)| < niax^gydistG(i),s)<6 -^(^)' that is 
we assign each vertex to an element of S within distance 5, but we increase the capacity of each s to the 
maximum capacity reachable within distance 6. In the second step we reassign vertices in such a way, that 
the maximum distance is at most 11, and at the same time capacity constraints are satisfied. 

Let T be any spanning tree of G^[S], rooted at an arbitrary vertex r £ S. Consider the following 
procedure, which constructs a function u : S ^ M_|_ U {0} assigning a tentative fractional number of centers 
to open in each of the vertices in S. 

1. Initially set u{s) = for each s £ S. 

2. For each vertex v £ V, if distG'(u, S) < 1, then by the fact that S is independent in G^ there exists 
exactly one vertex s^, G S, such that distG(f , s^) < 1, otherwise (when distG'(f , S) = 2), as s„ set 
any vertex from S within distance 2 from v. 

3. For each vertex v £ V increase u{sy) by (note that after this operation for each s G 5 we have 

u{s) > 1 and J2s&s'^(^) ~ ^)- 
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4. For each vertex s G 5 in a bottom-up order with respect to T, let x = u{s) — \u{s)'\ , decrease u{s) 
by X and increase u{p{s)) by x, where p{s) is the parent of s in the tree T (assume that p{r) = r). 

The last step of the above process ensures, that the function u has only integral positive values and moreover 
Sse5 ^(^) ~ ^- claim that if we open exactly u{s) centers in a vertex s, for each s G 5, then there 
exists an assignment (p : V ^ S, such that for each v ^ V we have distG(f , cl){v)) < 5 and for each s G 5 
we have < max^gyjig^^^j, g)<g L{v). Note that if we know that it exists, then we can find it in 

polynomial time by using maximum flow computation. To show the existence of (p, observe that it is enough 
to show a function f : V x S ^ U {0}, satisfying: 

• for each v e V we have J2seS /(^' ^) ~ Vv' 

• for each s e V we have Ylv&v fi"^^ ^) ~ "^i^)' 

• for each v e V,s £ S if f{v, s) > 0, then distG(v, s) < 5. 

Less formally, the function f{v, *) : 5 — )■ M+ U {0} is a distribution of the value y^, among vertices of S 
within distance 5. Such a function guarantees that we can fractionally cover all the vertices within distance 5 
when we open an integral number of centers u{s) in each vertex s £ S, therefore we can also cover vertices 
of V integrally, which proves the existence of the desired assignment (p. Observe, that we can construct the 
function / while performing the bottom- up process in Step|4j where for each vertex v £ V we split the value 
between f{v, s^) and f{v,p{sy)). It is always possible, because while going up the tree we can ensure 
that the part of y^, that was assigned to f{v,p{s.u)) remains in p{s.u), since we send up u{s) — \u{s)'] < 1 
and after Step |3] we have u{s') > 1 for each s' G S. Consequently, we know that there exists the desired 
assignment (p. 

Finally, we construct the final assignment (pQ : V ^ S' as follows. For each s G 5 let 

sl = arg max L{v) , 

v(^V,distc: (s ,v)<6 

open u{s) centers in and assign all the vertices of (p~^{s) to sl in 0o- By the properties of (p we infer that 
in this way we obtain a multiset S' of exactly k centers and an assignment (pQ, satisfying distG(f , <Po{v)) < 
11 for each v £ V and \(p~^{s')\ < L{s') for each s' £ S. □ 

Corollary 33. The integrality gap of LP 1 for connected graphs is at most 11 and there is a 11-approximation 
algorithm for connected graphs. 

Since Theorem|5]shows that a c-approximation algorithm for connected graphs implies a c-approximation 
algorithm for general graphs, the proof of Theorem|4]follows. 

4 Uniform capacities 

Here, we prove Theorem [8j i.e. we show 5 and 6 upper bounds for the integrality gap of LPl for uniform- 
soft-capacities and uniform-hard-capacities respectively, which is a counterposition to lower bounds of 4 
and 5 as stated in Theorem |7l 

The algorithm of KhuUer and Sussmann Il22l gives 6- and 5-approximation algorithm for the uniform 
capacitated A;-center problem, for hard and soft capacities respectively. It is possible to reformulate this 
algorithm to make it a rounding algorithm for LPl. However to avoid rewriting the whole algorithm, we 
show that if the algorithm of [22J does not produce a solution (which means that there is no solution which 
uses at most one hop), then it produces a witness, showing that there is no solution using k vertices and 
covering within distance of at most one hop. For a given instance {G, k, L) of uniform capacitated A;-center 
problem if the algorithm of KhuUer and Sussmann does not produce a solution, both in the case of hard and 
soft capacities, then it creates a set Vb ^ ^(C), such that: 
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(i) for any pair of vertices u,v £ Vq, the distance between u and t> in G is at least 3, 

(ii) \Vo\ + ^ > A:, where V = {v ^ V{G) : V„evbdistG(u, > 3}. 

Using the properties (i) and (ii) of Vq we prove that there is no feasible solution for LPl. By (i) and by 
constraints (|4]) and (|6]l of LPl we have: 

l^ol = XI X] ^u,v<'^ ^ yu= ^ Vv, (18) 

veVoueNolv] veVoueNo[v] veNc[Vo] 

where the last equality follows from (i). We lower bound the sum of y-values of vertices of A'^G'[y] = 
V\Ng[Vo] as follows. 

vgV ueNalv] w^Ng[V'] 

The first equality follows from constraint Q and ([6]) of LPl, whereas the inequality follows from Q of 
LPl. Therefore by ( [TSj ), ( 19 1 and (ii), we infer constraint ([T]l of LPl is violated. 

ly I 

v€V v£Ng[V'] v£Ng[Vo] 

Hence there is no feasible solution for LPl, which finishes the proof of Theorem [8] □ 



5 Integrality gap lower bounds 

In this section we present lower bounds on the integrality gap of LPl, i.e. prove Theorem[7] 

We start with a construction for uniform hard capacities. Let > 24 be an integer and set the uniform 
capacity as L = A; — 1. Let G' be a graph which consists of two adjacent vertices a, b together with L + 2 
vertices adjacent to both a and b. Let G be a graph composed of: 

• k — 6 copies of the graph G' denoted as for i = 1, . . . , A: — 6, 

• a single star Sk-e, rooted at r, which has exactly k — 6 leaves denoted as ^i, . . . , £k-6-> 

• k — 6 vertices for i = 1, . . . , A; — 6, 

• A; — 6 edges Xi,ii, 

• for each i = 1, . . . , k — 6 two edges between ai,bi G V{Gi) and Xj. 

Note that the graph G is connected. Observe that by setting the y-value equal to 1 for the vertex r, and 
Vai = Ubi = — ^ QSich i = 1, . . . , A: — 6, we can set x variables appropriately to obtain a feasible 
solution to the LP relaxation since the sum of y variables is equal to 1 + 2(A: — 6) ^jj^ = 1 + < 

Let us assume that there is a hard-capacitated A;-center in the graph G^. By Bi (for 1 < i < A; — 6) 
let us denote the set of vertices of Gj of degree two. Consider any 1 < i < A: — 6, and let us focus 
on how vertices of Bi can be covered. Observe that the only vertices in G that are within distance 4 
from Bi are V{Gi) U {xi,£i,r}. Since \Bi\ = L + 2 > L, at least two centers are opened in the set 
V{Gi) U {xi,ii,r}, which mean that at least one center is opened in V{Go) U {xi,£i,r}. Furthermore 
for different values of i, the sets V{Gi) U {xi,ii} are disjoint. Moreover if only one center is opened in 
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L + 2 



Figure 6: The graph G used in the proof of Theorem [7] 



the set V{Gi) U {xi,li}, then at least two vertices of Bi have to be covered by r. Hence even if a center 
is opened in r, then for at least (A; — 6) — ^ values of i we need to open at least two centers in the set 
V{Gi)yj {xi,li}. Moreover we need to open at least one center in each of the sets V{Gi)U {xi,li}, since 
otherwise L + 2 vertices can not be covered. In total there are at least (A; — 6) + ((fc — 6) — ^) centers 
opened, but (A: - 6) + ((A; - 6) - ^) = A + (A - 12) - ^ = A + > A since A > 24, a contradiction. 

For the uniform-soft-capacitated case the analysis is even simpler since if G^ admits a soft-capacitated 
A-center then without loss of generality there are at least two centers open in each of the sets V{Gi)VJ{xi,li\. 

To prove the lower bound for non-uniform capacities, in the above example we change capacities to zero 
for all the vertices except r and each of the 2(A — 6) vertices ai, hi. Note that those are the only vertices with 
non-zero y- value in the constructed feasible solution for LPl. Moreover we have shown that in any integral 
solution there exists a vertex of Bi, which has to be covered by a vertex of V{Gj) U {xj,lj}, for j / i, and 
since the only vertices of non-zero capacity in this set are the vertices Oj, hj, we infer, that some vertex of 
Bi has to be covered by a vertex at distance at least 7. Consequently the proof of Theorem[7]follows. □ 

6 (3 — e) -approximation hardness 

In this section we show, that it is not possible to approximate the A-center problem with non-uniform capac- 
ities in polynomial time with approximation ratio (3 — e) for constant e > 0. 

Theorem 34. If there exists a (3 — e)-approximation algorithm for the k-center problem with non-uniform 
capacities, where all the non-zero capacities, then P = NP. The theorem applies both to hard and soft 
capacitated version. 

Proof. Let us assume that a (3 — e) -approximation algorithm exists. We present a reduction from the EXACT 
Cover by 3-Sets problem, which is NP-complete. 

Exact Cover by 3-Sets 

Input: A set system (3", U), where each set in 3" has exactly 3 elements. 

Question: Does there exist a subset 3"' C 3", such that each element of the universe U belongs to exactly 
one set in 3"'. 

Let / = (J, U) be an instance of EXACT COVER BY 3-Sets. As a graph G we take the bipartite graph 
(JU5,^),whereS = [jfJ^^Ui,Ui = {ui : u £ U} and E = {Sm : 5 G J A 1 < i < | J| + 1 A -u G 5}. 
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That is G is an incidence graph of the set system (3", U), where the universe is repUcated I?"! + 1 times. 
Additionally for each 5 G 3" we add to the graph a vertex xs, which is adjacent to S, and has exactly 
3|3"| + 1 pendant neighbors (see Fig.|7]). We set a capacity function L : V{G) — Z+ U {0} as follows: 

• For each ui G Ui, L{ui) = 0. 

• For each 5 G J, L{S) = L{xs) = ^'3'\+ 3. 

• For each pendant vertex v, which is adjacent to xs, for G 3", we set L{v) = 0. 



3|:^| + 1 




L(v) = 

L{v) = 3|J^| +3 

L{v) = 3|7'| +3 



L{v) = 



U2 



Figure 7: The graph G constructed for the set system (3", [/), where 3" = {{A = {a,b,c},B 
{c, d, e}, G = {d, e, /}}} and U = {a, b, c, d, e, /}) 



Note that all the vertices, which have non-zero capacities, have exactly the same capacity. Let /' = 
(G, L, A; = |3"| + |[/|/3) be an instance of the capacitated A;-center problem. In what follows we prove: 

(i) If / is a YES-instance, then there exists a set Vq C V{G), such that | Vo| = k, together with a function 
4) : V ^ Vq, which satisfies V„gy(G)distG(w, < 1 and V^gvbl0~^('i^)l < -^(^)- Less formally, 
Vq is a solution for the capacitated A:-center problem, which opens at most one center in each vertex. 

(ii) Let Vq be a multiset containing exactly k vertices of V{G), such that there exists a function (j) '■ 
V ^ Vq, which satisfies y^^y(^Q^distG{v,4>{v)) < 3 and V^gvii < L{v). Then one can in 
polynomial time construct a set 3"' C 3", such that each element of the universe U belongs to exactly 
one set in 3"'. Intuitively, given a solution to the capacitated /c-center problem which uses distances at 
most two, and potentially opens more than one center in a vertex, one can in polynomial time construct 
a solution for the instance /. 

Observe that having (i) and (ii) suffices to prove the theorem, since if J is a YES-instance, we can obtain a 
solution to /, by constructing the graph G, running the (3 — e) -approximation algorithm on {G, L, k), which 
by (i) returns a set Vq, which by (ii) we transform to a solution for / in polynomial time. 

First we prove (i). Let 3^' be a solution for the instance I. We take Vq = 3"' U {xs : 5" G 3"}, that is 
\U\/3 sets from 3", and all the vertices xs- Note that |Vo| = \U\/3 + |3"| and each vertex in Vq has capacity 
exactly 3|3"| + 3. Observe that if each vertex S" G 3~' covers each of the |3"| + 1 copies of the three elements 
in S, while each vertex xs covers itself, all the 313"! + 1 pendant vertices, and the vertex S, then we obtain 
the desired function (p. 
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Next we prove (ii). Observe, that \V{G)\ = (|J| + 1)\U\ + |:r|(3|J| + 3) = (3|J| + 3)(|[/|/3 + |^|), 
and hence in the multiset Vq there are only vertices of capacity 3|3"| + 3, since otherwise the total capacity 
of the vertices in Vq would be smaller than |y(G)|. Moreover, without loss of generality we can assume 
that the multiset Vq contains only vertices 5 G 3", since we can always replace a vertex X5 by S without 
exceeding distance two, in the function 0. Furthermore pendant vertices of any xs are covered by vertices 
at distance at most two, that is by the vertex S. Therefore each 5 E 3" appears in the multiset Vq at least 
once. Let 3"' C 3" contain exactly those sets 5 G 3", for which the vertex S belongs to Vq more than once. 
Note that since |Vo| = |3"| + |f/|/3, the set 3~' contains at most \U\/3 sets. Consequently to prove (ii) it is 
enough to show that each element of the universe U belongs to at least one set in 3"'. 

Let Vq C Vq be the set containing the first copy of each G 3" in Vo» namely Vq = 3". Without loss of 
generality, for each S G Vq the set (j)~^{S) contains xs and all its pendant neighbors, since none of those 
vertices can be covered by S' / S, because the distance would be at least three. Consequently for each 
S e Fq' we have I^^^S) n (U2i^^ Ui)\ < L{S) - (3|:r| + 2) = L Since l^o'l < | J| < |3"| + 1, there 
exists an index I < io < + 1, such that for each S G Vq we have (f)^^{S) n C/j = 0. Therefore each 
vertex Uj G Ui is covered by some set S" G 3"' that it belongs to, which proves (ii) and finishes the proof of 
Theorem [34l 

□ 

7 Conclusions and open problems 

We have obtained the first constant approximation ratio for the fc-center problem with non-uniform hard 
capacities. The approximation ratio we obtain is in the order of hundreds (however we do not calculate 
it explicitly), so the natural open problem is to give an algorithm with a reasonable approximation ratio. 
Moreover, we have shown that the integrality gap of the standard LP formulation for connected graphs in 
the uniform capacities case is either 5 or 6, which we think might be an evidence, that it should be possible 
to narrow the gap between the known lower bound of (2 — eps) and upper bound 6 in the uniform capacities 
case. 
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